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Abstract 

We estimate the decay constants f n and fx of the light mesons within a 
relativistic quantum-field model of interacting quarks and gluons confined an- 
alytically. The necessary physical parameters, the quark masses m u and m s , 
the coupling constant a s and the confinement scale A have been obtained from 
our previous investigation on the meson ground states, orbital and radial excita- 
tions by using the ladder Bethe-Salpeter equation. Our model provides a solid 
framework to compute the meson spectra, the lowest glueball state and decay 
constants of light mesons from the basic principles of QCD and QFT. 
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^ ! 1 Introduction 



The pseudoscalar mesons, especially the n and K, are studied extensively to under- 
O 1 stand the under-structures from nonperturbative QCD (e.g., [1]). The decay constants 
of the light mesons are one of the basic parameters in the particle physics and form an 
Ph! extremely important test ground for any nonperturbative QCD methods. Particularly, 
the pseudoscalar-meson decay constants fp play an important role in hadron phe- 
nomenology (e.g., the leptonic decay width is proportional to fp and the V us element 
is extracted from the ratio fx/ fn) an d can be used as a test ground for models. 

A number of theoretical approaches obeying an accuracy comparable to the uncer- 
tainties in the experimental data, has relied on models containing too many, in some 
cases even nonphysical, parameters. On the other hand, the scaling property of QCD 
implies that all characteristics of hadrons must depend on a few global parameters. 

Among the significant theoretical approaches, the Bethe-Salpeter equation (BSE) is 
a conventional approach in dealing with the two-body relativistic bound state problems 
[2] and its solutions give useful information about the under-structure of the mesons. 
Many analytical and numerical calculations indicate that the ladder BSE with phe- 
nomenological models can give model independent results successful descriptions of 
the long distance properties of the low energy QCD and the QCD vacuum [3, 4, 1]. 

Below we adopt a relativistic quantum field model of interacting quarks and gluons 
under the analytic confinement, solve the BSE in the one-gluon exchange approximation 
for the quark-antiquark bound states and estimate the decay constants f n and fx- The 
model relies on only necessary physical parameters {mf,a s ,A} which have been fixed 
to fit the meson masses in the ground states, orbital and radial excitations as well as 
the lowest-state glueball mass [5, 6]. 
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2 Two-quark Bound States 



In the conventional quark model mesons are (qq') bound states of quarks and antiquarks 
(the flavors may be different). We deal with the pseudoscalars P (J PC = CT + ) and 
vectors V (1 ), the most established sectors of hadron spectroscopy. 

Consider a Yukawa-type interaction of quarks and gluons with the Lagrangian 

c = (vs- 1 *) + \ (0 d-v) + g (K^%%<p;) , (i) 

where ^ a a f(x) is the quark field, <f>^(x) is the gluon field, g is the coupling strength, t a 
are the Gell-Mann matrices and {a, a, /} are the color, spin and flavor indices. 

The conventional QCD encounters a difficulty by defining the explicit quark and 
gluon propagator at the confinement scale [7]. Particularly, the infrared form of the 
gluon propagator is still a controversial aspect [8], there are various versions, with 
different behaviors for the infrared region: finite, zero and more divergent that 1/k 2 , 
each one with its advantages and inconveniences [9, 10, 11, 12]. 

Nevertheless, the matrix elements of hadron processes at large distance are in fact 
integrated characteristics of the quark and gluon functions and taking into account 
the correct global symmetry properties and their breaking by introducing additional 
physical parameters may be more important than the working out in detail (e.g., [13, 
4, 14]). Besides, the bound states of quarks and gluons may be found as the solution 
of the BSE in a variational-integral form [5, 6] that is low sensible on tiny details of 
integrands. 

We consider the simplest effective quark and gluon propagators: 

$iw^ Kl] {w + m [ 1 + l^M]} a( s f p 2 + m 2 \ 
W= S' -> eXp {-^A^| ' 

DZ(x) = ^ exp j-^} = D(x) , (2) 

where p = p^^, m - the quark mass and < iv(m) < 1 [6]. These entire analytic 
functions in the Euclidean metric keep the spirit of the analytic confinement [15, 16, 4] 
and imply that each isolated quark and gluon is confined in the background gluon field. 
The meson mass is derived from [6]: 

l + A A r(A4) = 0. (3) 

where the two-quark bound states are obtained by diagonalizing the quadratic part in 
the meson partition function on the full orthonormal system {U^}. This is nothing 
else but the solution of the appropriate BSE (in the one-gluon approximation): 

jj dxdy U N {x) {l + g 2 jD(x) Tljj,(p 2 , x, y)^D{^ U N ,{y) = 5 NM '[l + \ N {-p 2 )] , (4) 
where, the two-point function is introduced: 

UjAp 2 , x,y)=f e~ ik ^Tr [TjS (k + ^p) Tj,S (k - ^p) 

and Tj = {ry 5 , ry M } for J = {P, V}. 
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We have solved Eq.(3) within a variational-integral approach and have found that 
the following values of free parameters: 



a. 



m. 



= 0.186, 
188 MeV . 



A = 730 MeV 



m r 



646 MeV . 



m {u, d} 

m b -- 



-- 170 MeV, 
4221 MeV. 



(5) 



are optimal for the meson ground states, orbital and radial excitations [6]. Our esti- 
mates for the P and V meson masses in the ground state (Fig.l) were in good agreement 
with the experimental data [17] and have a small experimental error ( 1 -f- 5 per cent) 
in a wide range ~ 140 4- 9500 MeV. 
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Figure 1: Estimated meson masses (dots) compared with experimental data (dashes) 
The lowest-state (scalar) glueball mass has been derived analytically [6]: 

2vr(2 + v^) 2 



Ml 



2A 2 In 



a 



upp 



O.s 



(1745 MeV) z 



a 



Upp 



27 



(6) 



Our estimate (6) is close to the result Mq = 1710 ± 100 MeV due to the quantized 
knot soliton model [18] and quenched lattice estimate [19] and in agreement with the 
predictions expecting non-gg scalar object in the range ~ 1600 4 1800 MeV [20, 21]. 

The estimated value a s = 0.186 coincides with the latest experimental data ot s t expT \ ^ 
0.10 -=-0.35 [17] and the prediction of the quenched theory ot s ^ quenc h) ~ 0.195 [22]. Note, 
this relatively weak coupling justifies the use of the one-gluon exchange mode in our 
consideration. By analysing the glueball mass formula (6) one obtains the restriction 
®s < Oi upp ~ 3.24124 that coincides with the two-loop value of the freezing coupling 
constant [23]. 

A rough estimate of the hadronization distance 



Jd 4 xx 2 D(x) 8 / 1 



hadr 



Jd 4 xxD(x) A 2 V250MeV 



shows that the confinement radius r con f ~ 1/ A and r^adr are comparable values. 
Below, we extend our consideration to the light meson decay constants. 
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3 Light Meson Decay Constants 

The decay constant fp is denned by the matrix element: 

ifpp ll = <0|J M (0)|E/(p)) , 



(7) 



where J M is the axial vector current and U(p) is the normalized state vector. 

There have been made many attempts to extract fp from meson spectroscopy. 
Particularly, the decay constant f n + for 7r + -meson is determined from the combined 
rate for tt + — > + (j, + v l /y. With the recent experimental data one obtains [17]: 

f* xp = 130.7±0.1±0.36MeV, f e K xp = 159.8±1.4±0.44MeV , fT 1 17 V ^ 1.22±0.02 . 

Recent lattice QCD model with exact chiral symmetry estimates [24]: 



latt 
K 



152±6± lOMeV. 



The unquenched lattice QCD calculation predicts [25] 



7r Junquen 



1.198 ±0.019. 



Following our approach, we estimate fp as follows: 
g f dk 



ifp Vy. 



dx e- ikx U(x) ^D{x)Ti [i^S (k + f^p) i lsl ^S (k - fi 2 p) } . 



6 J (2?r) 4 

It is reasonable to choose a normalized trial function [6] : 
U{x) ~ y/D(x) exp {-aA 2 x 2 /4:} , f dx \U(x)f 



a > . 



(8) 
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Figure 2: Evolution of f n and fx with A. Dashed lines depict the experimental data. 
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Omitting details of intermediate calculations write the final expression 



Rx s (mi + m 2 ) 



127r 3 / 2 (mi/A) 2 (m 2 /A)< 



exp 



ml + m\ 
2A 2 



Ml 



max < (3c — 1)(1 

1/2<c<1 1 



(m 1 + m 2j 

1m\m<i + c(mi — m 2 ) 2 
c) ; — exp 



(mi + m 2 ) 



- 1 



cMl 



(9) 



mi — m 2 



4 A 2 \mi + m 2 



The solution of Eq.(9) for different values of A is plotted in Fig. 2 in comparison 
with the experimental data [17]. In the region 650 -j- 760MeV we underestimate f% xp 
and overestimate f^ v ■ Evolutions of fx/ fn and a s with A are depicted in Fig. 3. 




Figure 3: The ratio fx/ fn and the coupling strength a s versa the confinement scale A. 
For consistency, we use the parameters (5) and obtain: 

f w = 118MeV , f K = 176MeV . 
Our estimates lie near to the latest data [17, 24, 25]. 

In conclusion, we have considered a relativistic quantum field model of interacting 
quarks and gluons under the analytic confinement and solve the Bethe-Salpeter equa- 
tion in the one-gluon exchange approximation for the quark- ant iquark (qq') bound 
states. Despite the simplicity, this model resulted in a quite reasonable sight to the 
underlying physical principles of the hadronization mechanism and gave an accurate 
description of the ground states, orbital (£ > 0) and radial {n r > 0) excitations of 
mesons in a wide range of mass (up to 9.5 GeV) and acceptable values of the light 
meson decay constants {/,,-, fx}- Our model provides a solid framework to compute 
the meson spectra, the lowest glueball state and decay constants of light mesons from 
the basic principles of QCD and QFT. 
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